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Optica]  signal  transmission  is  an  important  component  of  nu¬ 
merous  underwater  applications,  including  visibility  and  electro- 
optical  (EO)  communication.  In  addition  to  the  well-studied  effect 
of  particle  backscatter,  underwater  optical  signal  transmission  can 
be  limited  by  changes  in  the  index  of  refraction  (IOR)  due  to  small- 
scale  variations  in  temperature  and  salinity,  sometimes  called  "op¬ 

tical  turbulence”.  These  variations  in  IOR,  which  are  associated 
with  oceanic  turbulence,  can  lead  to  the  blurring  of  an  underwater 
optical  target,  particularly  at  high  spatial  frequencies,  thus  reduc¬ 
ing  target  detail.  The  201  1  BahamasOpticalTurbulence  Experiment 
(BOTEX)  was  conducted  to  investigate  this  impact  of  turbulence 
on  underwater  optical  signal  transmission.  Investigating  naturally 
occurring  “optical  turbulence”  requires  a  platform  held  at  depth, 
capable  of  concurrent  measurements  of  optical  impairment  by  tur¬ 
bulence,  which  requires  a  significant  optical  path  length,  as  well 
as  associated  physical  and  optical  background  conditions  of  the 
ambient  environment.  Our  novel  platform  consisted  of  a  high¬ 
speed  camera  and  optical  target  mounted  on  a  5m-long  frame, 
along  with  several  Nortek  Vector  Acoustic  Doppler  Velocimeter 
(ADV)  and  PME  Conductivity-Temperature  (CT)  probes,  to  estimate 
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turbulent  kinetic  energy  and  temperature  variance  dissipation 
rates  experienced  by  the  frame.  Data  on  the  background  turbu¬ 
lence  was  collected  with  a  Rockland  Oceanographic  Vertical  Mi¬ 
crostructure  Profiler,  to  aid  in  analysis  and  guide  error  estimates  of 
the  ADV/CT  measurements.  This  study  was  the  first  effort  attempt¬ 
ing  to  collect  turbulence  measurements  on  a  frame  designed  for 
the  investigation  of  the  effect  of  density  microstructure  variations 
on  optical  signal  transmission  in  the  open  ocean.  Our  results  high¬ 
light  the  numerous  challenges  associated  with  studying  this  phe¬ 
nomenon  in  the  dynamic  oceanic  environment.  Here,  we  present 
the  interpretation  of  the  high-resolution  velocity  and  temperature 
measurements  collected  on  the  frame  and  discuss  the  associated 
difficulties.  Despite  the  numerous  challenges,  the  investigation  of 
the  effect  of  microstructure  on  underwater  optics  is  needed  for 
efforts  aimed  at  mitigating  the  impact  of  “optical  turbulence”  on 
underwater  EO  signal  transmission  and  may  help  advance  optical 
methods  to  quantify  oceanic  microstructure. 

Published  by  Elsevier  B.V. 


1.  Introduction 

Optical  properties  of  coastal  and  open  ocean  water  and  their  impact  on  electro-optical  (EO)  im¬ 
agery  are  important  for  a  wide  range  of  applications.  Diver  visibility,  mine  detection,  search  and  res¬ 
cue,  and  optical  communications  all  depend  on  underwater  optical  signal  transmissions.  Underwater 
imaging  has  been  shown  to  be  affected  by  absorption  and  scattering  due  to  suspended  particles  (Hou 
et  al.,  2007).  However,  recent  evidence  suggests  that,  in  addition  to  particle  scattering,  underwater 
EO  imagery  may  also  be  significantly  affected  by  turbulent  temperature  and  salinity  microstructure 
in  the  water.  Turbulence  is  ubiquitous  in  the  world’s  oceans  and  fresh-water  environments  and  the 
associated  small-scale  temperature  and  salinity  variations  can  lead  to  localized  changes  in  the  index 
of  refraction,  which  in  turn  impacts  underwater  optical  properties.  It  has  been  demonstrated  that 
this  so-called  “optical  turbulence”  can  be  a  limiting  factor  in  natural  environments,  affecting  opti¬ 
cal  signal  transmission  that  impact  various  applications,  from  diver  visibility  to  active  and  passive 
remote  sensing  (Hou  et  al.,  2012b;  Bogucki  et  al.,  2004).  The  recent  Skaneateles  Optical  Turbulence 
Experiment  (SOTEX)  conducted  in  a  lacustrine  environment  confirmed  that  optical  turbulence  may 
significantly  contribute  to  image  degradation  (Hou  et  al.,  2012b;  Woods  et  al.,  2011).  To  quantify  this 
impact  of  turbulence  on  optical  signal  transmission  in  an  oceanic  environment,  the  Bahamas  Optical 
Turbulence  Experiment  (BOTEX)  was  conducted  in  the  summer  of  2011  (Fig.  1)  (Hou  et  al.,  2012a). 
Measuring  turbulence  in  the  ocean  has  long  been  known  to  be  a  challenging  task,  and  thus  study¬ 
ing  naturally-occurring  optical  turbulence  in  the  open  ocean  is  rife  with  difficulties.  Among  these  are 
the  deployment  strategy,  as  well  as  ambient  atmospheric  and  oceanographic  conditions  that  can  also 
affect  the  measurement  method.  The  goal  of  this  paper  is  to  describe  and  address  the  challenges  en¬ 
countered  while  working  to  measure  turbulence  parameters  with  an  experimental  array  that  is  geared 
towards  optics  measurements  and  different  from  the  traditional  turbulence  profiling  or  fixed  mooring 
approach.  Since  turbulence  is  patchy  and  highly  variable  in  both  space  and  time,  it  is  important  to  have 
a  type  of  turbulence  sensor  on  the  same  structure  that  is  also  measuring  information  on  the  optics. 
This  limits  what  turbulence  sensors  can  be  used.  Shear  probes,  such  as  used  in  traditional  Vertical  Mi¬ 
crostructure  Profilers  (VMP)  and  generally  deployed  at  a  free-fall  from  the  vessel,  would  not  be  useful 
on  a  tethered  frame  drifting  with  the  boat.  The  requirements  on  optical  path  length  further  complicate 
the  setup.  An  optical  path  length  of  several  meters  is  desirable  in  particle-laden  natural  waters,  since 
the  effect  of  image  distortion  due  to  turbulence  is  cumulative  and  often  less  pronounced  than  that 
of  blurring  due  to  particle  scattering  (Duntley,  1963;  Preisendorfer,  1986;  Hou  et  al.,  2007).  This  puts 
certain  size  requirements  on  the  frame  to  measure  the  optics,  which  in  turn  limits  the  ability  to  let 
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Fig.  1.  Location  of  stations  during  BOTEX. 

the  array  free-fall  through  the  water  column.  In  the  present  study,  we  settled  on  the  use  of  Acoustic 
Doppler  Velocimeters  (ADVs)  and  high-resolution  thermistors  (PME  Conductivity-Temperature  (CT) 
probe)  affixed  to  the  frame.  However,  the  use  of  ADVs  to  infer  turbulent  kinetic  energy  dissipation 
rates  (TKED)  is  known  to  be  inherently  more  successful  in  more  energetic  environments.  In  more  qui¬ 
escent  waters,  the  signal  to  noise  ratio  of  the  sensor  may  reach  its  lower  limit.  Since  noise  from  plat¬ 
form  motion  and  vibrations  can  complicate  the  estimation  of  turbulence  properties  from  the  ADV/CT 
package,  complementary  measurements  of  ambient  turbulence  levels  with  a  dedicated  microstruc¬ 
ture  profiler,  such  as  the  VMP,  are  an  invaluable  tool  to  aid  in  data  analysis  and  guide  error  estimates. 

This  paper  describes  our  instrument  setup,  the  procedures  used  to  calculate  turbulence  parameters 
from  the  measurements,  as  well  as  the  challenges  encountered  with  our  system.  We  also  briefly 
discuss  the  concurrent  optical  observations  and  their  relation  to  the  observed  turbulence  parameters. 
A  more  detailed  analysis  of  the  optical  measurements  taken  during  this  field  experiment  will  be  the 
subject  of  future  work. 

2.  Methods 

The  goal  for  the  measurements  obtained  during  this  study  was  to  be  able  to  image  across  a  signif¬ 
icant  optical  path  length  while  at  the  same  time  collecting  data  on  the  turbulence  properties  of  the 
ambient  water.  Deploying  a  free-falling  platform  similar  to  a  traditional  turbulence  profiler  would 
severely  limit  the  optical  path  length  and  hinder  the  ability  to  image  across  naturally  occurring  opti¬ 
cal  turbulence,  such  as  in  the  thermocline.  Imaging  across  such  naturally  occurring  density  variations 
in  index  of  refraction  (IOR)  is,  of  course,  also  complicated  if  the  frame  is  tethered  to  a  drifting  vessel 
and  affected  by  cable  motion  and  vibrations,  where  the  frame  itself  may  induce  turbulence.  The  ap¬ 
proach  chosen  was  a  compromise,  providing  a  significant  optical  path  length  along  the  measurement 
structure  -  a  5  m  long,  rigid  aluminum  frame,  termed  the  Image  Measurement  Assembly  for  Subsur¬ 
face  Turbulence  (iMast),  -  while  simultaneously  collecting  high-resolution  velocity  and  temperature 
data  for  estimation  of  turbulent  properties  of  the  water.  A  similar  frame  was  used  during  the  same 
field  exercise  to  study  laser  sensing  of  mixed  layer  turbulence  (Dalgleish  et  al.,  2013),  however,  these 
measurements  did  not  collect  any  turbulence  data  directly  on  the  frame,  but  rather  compared  the 
optics  results  to  the  background  values  from  the  VMP.  Depending  on  sea  state  and  atmospheric  con¬ 
ditions,  as  well  as  due  to  the  patchy  and  time-variable  nature  of  turbulence,  such  an  approach  omits  an 
important  piece  of  information,  the  turbulence  encountered  at  the  time  of  the  optical  measurements, 
as  well  as  frame-induced  turbulence. 

Thus,  in  order  to  characterize  the  velocity  and  temperature  microstructure  experienced  by  the 
optics  on  our  frame,  high-resolution  velocity  and  temperature  measurements  were  obtained  with 
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Fig.  2.  Schematic  of  iMast  deployment  during  BOTEX  (top)  and  photograph  of  iMast  on  deck  showing  the  ADV/CT  probes 
mounted  on  the  frame,  the  high-speed  camera  and  the  iPad  as  optical  target  (bottom). 


a  Nortek  Vector  Acoustic  Doppler  Velocimeter  (ADV)  and  PME  Conductivity-Temperature  (CT)  mi¬ 
crostructure  probe  simultaneously  with  measurements  of  optical  target  clarity,  i.e.,  a  high-speed  cam¬ 
era  was  continuously  imaging  an  optical  target  across  a  region  of  turbulence  or  quiescent  water.  Since 
turbulent  properties  from  the  ADV/CT  data  are  difficult  to  estimate  and  are  sensitive  to  contamina¬ 
tion  by  platform  and  cable  motion,  background  turbulence  levels  were  also  assessed  with  a  traditional 
turbulence  profiling  instrument.  This  Vertical  Microstructure  Profiler  (VMP)  measures  velocity  shear 
and  temperature  while  freely  falling  through  the  water  column  (VMP  is  manufactured  by  Rockland 
Oceanographic  Services  Inc.,  Canada  and  similar  to  a  profiler  described  in  Wolk  et  al.,  2002).  The  TKED 
and  temperature  variance  dissipation  rates  (TD)  were  calculated  from  the  ADV/CT  measurements  via 
spectral  fitting  to  Kolmogorov  spectra  (for  velocity)  and  Batchelor  spectra  (for  temperature)  and  com¬ 
pared  to  TKED  and  TD  estimates  from  the  VMP  (Oakey,  1982). 

2.1.  Data  collection 

The  ADV  and  CT  probes,  along  with  a  high-speed  imaging  camera  on  one  end  and  optical  tar¬ 
gets  (active  and  passive)  on  the  other,  were  mounted  on  the  rigid  aluminum  frame  of  the  iMast.  The 
iMast  structure  was  deployed  from  the  R/V  Walton  Smith’s  A-frame  and  lowered  to  a  set  of  predeter¬ 
mined  depths  (Fig.  2)  (Hou  et  al.,  2012b).  The  frame,  tethered  to  the  ship,  was  held  at  these  depths 
for  10-15  min  long  “depth  pauses”,  to  collect  time  series  of  high-resolution  velocity  and  tempera¬ 
ture/conductivity  for  the  subsequent  estimation  of  TKED  and  TD  rates.  Data  were  collected  at  a  sam¬ 
pling  frequency  of  32  Hz.  The  pauses  were  taken  at  depths  where  optical  turbulence  associated  with 
strong  temperature  gradients  was  expected  to  be  most  pronounced,  i.e.  at  the  bottom  of  the  mixed 
layer  and  just  below  the  thermocline,  although  results  from  the  study  using  a  longer  frame  during  the 
same  field  experiment  showed  significant  cumulative  effects  across  most  of  the  upper  layer  (Dalgleish 
et  al.,  2013).  For  comparison,  additional  data  were  collected  at  depths  exhibiting  low  temperature 
variability.  The  estimation  of  TKED  and  TD  from  this  tethered  structure  is  complicated  by  platform 
movement,  and  hence  additional  data  on  the  background  turbulence  was  collected  with  the  VMP  ei¬ 
ther  immediately  preceding  or  succeeding  the  iMast  deployment.  These  measurements  to  describe 
the  ambient  turbulence  levels  were  taken  with  a  Rockland  VMP,  which  was  deployed  as  a  free-falling 
profiler  loosely  tethered  to  the  ship.  In  this  paper,  we  focus  on  the  data  collected  with  the  iMast; 
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Fig.  3.  Pressure  (left)  and  velocity  (right)  time  series  from  the  Florida  coast  (Station  1 ,  top)  and  the  Bahamas  (Station  3d,  bottom) 
illustrate  platform  motion  felt  by  the  ADV/CT  probes  on  the  iMast.  The  extent  of  vertical  movement  by  the  frame  is  similar  at 
both  stations,  but  the  period  of  motion  is  shorter  and  the  vertical  velocities  are  higher  at  the  Bahamas  station.  Meteorological 
conditions  were  recorded  as  calm  at  Station  1  (yet  with  strong  currents  and  drift),  and  as  rougher  with  pronounced  swell  at 
Station  3. 


data  collection  with  the  VMP  and  processing  of  its  data  follow  established  protocols,  more  details 
can  be  found  in  Oakey  (1982)  and  Hou  et  al.  (2012b).  In  addition  to  the  measurements  of  turbulent 
properties,  hydrographic  and  optical  properties  of  the  water  were  measured  with  a  SeaBird  Conduc- 
tivity-Temperature-Depth  (CTD)  and  a  WetLabs  ac-9  transmissometer  mounted  on  a  small  wire  cage 
and  deployed  from  the  vessel  at  each  station  before  or  after  the  iMast  or  VMP  deployments. 


2.2.  Calculation  of  TKED  and  TD 


2.2.1.  Turbulent  kinetic  energy  dissipation  rate— TKED 

The  estimation  of  TKED  rates  from  velocity  time  series  collected  with  ADVs  is  based  on  the  well- 
established  Inertial  Subrange  Dissipation  Method  (IDM)  (see,  f.  ex.,  Bluteau  et  al.,  2011)  and  is  based 
on  fitting  observed  velocity  spectra  to  the  theoretical  energy  spectrum  of  turbulence  described  by 
Kolmogorov  (see,  f.ex.,  Lesieur,  2008),  which  exhibits  a  —  |  slope  in  the  inertial  subrange  (ISR).  The 
spectral  fitting  requires  that  an  ISR  is  present  and  well-resolved  in  the  spectra.  Deployment  strategies 
thus  have  to  take  into  account  platform  motion  and  instrument  configuration,  as  these  can  introduce 
motion  that  may  cause  contamination  of  the  ISR  of  the  spectrum  (Fig.  3).  Resolving  the  ISR  also 
puts  a  minimum  requirement  on  the  record  length.  Ideally,  the  full  10-15  min  time  series  from  the 
depth  pauses  would  be  used  to  calculate  TKED.  However,  excessive  noise  in  the  data  due  to  platform 
vibrations  as  well  as  directional  changes  in  mean  flow  U  often  reduced  the  length  of  usable  data 
segments. 

Spectra  were  calculated  from  velocity  time  series,  and  Taylor’s  Frozen  Turbulence  hypothesis  (see, 
f.  ex.,  Tennekes  and  Lumley,  1972),  /<  =  is  applied  to  convert  the  frequency  spectra  to  wavenum¬ 
ber  spectra.  Here,  k  is  the  wavenumber  (in  rad  m-1 ),/  is  the  frequency  (in  s_1 ),  and  U  is  the  mean  flow 
past  the  sensor  (in  m  s-1).  A  wavenumber  range  representative  of  the  spectrum  in  the  ISR  is  chosen 
and  the  TKED  rate  e  is  estimated  from 

Eu(ki)  =  Cehi~5\  (1) 

where 

oo 

Eu(kj)  dkj  =  uf.  (2) 

Here,  E„  is  the  one-dimensional  velocity  spectrum  of  component  i  {x\  =  x  is  the  streamwise 
component,  x2  =  y  is  the  transverse  component,  and  x3  =  z  is  the  vertical  component,  in  a 
Cartesian  coordinate  system),  and  uf  is  the  variance  of  the  signal.  C  is  a  constant  equal  to  a 
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Fig.  4.  Energy  spectra  from  the  streamwise  velocity  component  at  Station  3  in  the  Bahamas  (Tongue  of  the  Ocean).  Note  the 
peak  at  wavenumber  k  %  10-20  rad/s  (corresponding  to  a  period  of  about  3-5  s)  in  the  unfiltered  spectra  (left),  which  is  due  to 
ship/platform  motion.  For  the  spectrum  shown  on  the  left,  U  was  chosen  to  be  the  larger  of  the  horizontal  velocity  components. 
The  filtered  spectrum  is  shown  on  the  right.  Here,  the  velocities  were  rotated  into  a  new  frame  of  reference  with  U  =  ( U , ,  0. 0). 
The  slope  in  the  inertial  subrange  of  the  spectra  is  in  agreement  with  the  Kolmogorov  -|  law. 

in  the  streamwise  direction  (x  and  ki)  and  to  |  ^|or  in  the  transverse  or  vertical  direction  (y,  z); 
or  ^  1.5  is  the  Kolmogorov  constant  (e.g.,  Rusello  and  Cowen,  2011).  Since  a  considerable  source 
of  uncertainty  in  the  estimation  of  €  is  introduced  through  the  choice  of  the  mean  flow  U  for  the 
application  of  Taylor’s  Frozen  Turbulence  hypothesis,  we  follow  Bluteau  et  al.  (2011)  and  rotate  our 
observed  velocities  (which  were  collected  in  XYZ -coordinates,  the  coordinate  system  of  thesensor, 
then  rotated  into  geophysical  coordinates)  into  a  new  frame  of  reference  where  U  =  (U i,  0,  0) 
(Bluteau  et  al.,  201 1).  If]  is  then  used  in  the  conversion  from  frequency  to  wavenumber  space.  For 
our  data,  the  wavenumber  range  to  fit  to  the  theoretical  spectrum  was  chosen  manually,  we  found 
this  approach  to  yield  a  more  robust  result  than  automating  the  routine,  especially  in  the  case  of 
noisy  spectra  (Fig.  4).  The  fit  to  the  data  was  accomplished  by  using  the  approach  introduced  by 
Bluteau  et  al.  (2011),  based  on  a  maximum  likelihood  estimator  (MLE)  technique,  previously  used 
for  Batchelor  spectral  fitting  (Ruddick  et  al.,  2000).  This  approach  also  allows  the  estimation  of  errors 
on  TKED  through  the  MLE  method  as  described  in  Bluteau  et  al.  (2011)  and  Ruddick  et  al.  (2000). 
These  errors  represent  the  quality  of  the  fit  of  the  observed  spectrum  to  the  theoretical  —  |  slope. 
They  do  not,  however,  adequately  describe  all  uncertainties  present  in  the  data  due  the  assumptions 
underlying  the  Frozen  Turbulence  Hypothesis  and  ISR  method. 

One  issue  that  arose  in  our  data  was  the  occurrence  of  so-called  “phase-wrapping”,  which  can 
occur  if  the  ambient  velocity  exceeds  the  maximum  velocity  that  can  be  measured  by  the  instrument, 
the  so-called  ambiguity  velocity.  This  maximum  measurable  velocity  of  the  ADV  is  determined 
by  an  instrument  setting  controlling  the  pulse  separation,  the  nominal  velocity  range,  and  applies 
to  the  beam  velocities  (Fig.  5).  Generally,  it  should  be  ensured  that  the  ambiguity  velocity  is  not 
exceeded  during  deployment,  however,  this  increase  in  nominal  velocity  range,  i.e.,  a  decrease  in  pulse 
separation,  also  leads  to  a  decrease  in  instrument  resolution.  Since  this  limit  on  measurable  velocity 
applies  to  the  beam  velocities,  phase-wrapping  can  occur  even  if  the  observed  earth-referenced 
velocities  do  not  exceed  the  estimated  maximum  velocity  range  (Rusello,  2009).  Data  with  mild  to 
moderate  amounts  of  phase-wrapping  can  be  corrected  for  in  post-processing.  The  correction  has 
to  be  applied  for  the  beam  velocities  and  requires  a  coordinate  transformation  if  data  is  collected 
in  XYZ  or  ENU  coordinates  (i.e.,  the  coordinate  system  relative  to  the  probe  head  or  the  geophysical 
coordinate  system,  respectively)  (Rusello,  2009).  Extreme  cases  of  phase-wrapping  can  be  too  noisy  to 
unwrap  successfully  and  may  introduce  spuriously  high  values  for  TKED  dissipation;  such  data  were 
rejected  in  our  analysis. 


2.2.2.  Temperature  variance  dissipation  rate—TD 

While  we  did  collect  conductivity  data,  we  chose  here  to  report  only  on  the  temperature  data.  It 
has  previously  been  reported  that  variations  in  the  index  of  refraction  in  the  ocean  are  mainly  driven 
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Fig.  5.  Phase-wrapping  occurs  when  the  so-called  beam  ambiguity  velocity,  uamb,  is  exceeded  (for  \uamb\  >  0.23  m/s  in  our 
case).  This  leads  to  spurious  velocities  in  the  probe-  or  earth-referenced  data  (left  panel,  blue  line).  The  ambiguity  velocity  is 
determined  by  an  instrument  setting,  which  controls  the  pulse  separation  in  pulse-coherent  Doppler  velocimeters.  When  data 
is  phase-wrapped,  the  correlation  score  between  successive  pulses  (should  be  >70%)  drops  (bottom  right  panel).  Such  data  has 
to  be  rejected  unless  it  can  be  corrected,  i.e.  “unwrapped”.  The  red  line  in  the  top  right  panel  shows  the  “unwrapped”  velocity, 
the  blue  is  the  phase-wrapped  data. 


by  temperature  fluctuations  (Bogucki  et  al.,  1998,  2004).  Furthermore,  we  felt  that  our  conductivity 
data  may  not  be  reliable  enough  to  properly  estimate  salinity  variance  dissipation  rates,  as  described 
in  Nash  and  Mourn  (1999),  especially  given  our  relatively  low  sampling  frequency. 

For  the  estimation  of  TD  rates  from  the  temperature  time  series,  we  followed  Mourn  and  Nash 
(2009),  fitting  temperature  gradient  spectra  to  the  theoretical  Batchelor  spectrum.  This  technique 
involves  fitting  the  temperature  gradient  spectra  in  the  dissipative  high-frequency  range  of  the 
spectrum.  To  obtain  temperature  gradient  spectra  from  the  temperature  time  series,  we  again  invoke 
Taylor’s  frozen  turbulence  hypothesis  and  use 


dT  _  1  dT 
~dx  ~  U~dt 


(3) 


and  k  =  f  /U  (note  that  k  has  units  of  cpm).  For  isotropic  turbulence,  the  dissipation  rate  of  the 
temperature  variance  can  be  described  as 


Xt  =  6DT 


# tx  (k)  dk, 


(4) 


where  DT  is  the  thermal  diffusivity,  and  &Tx  (k)  is  the  wavenumber  spectrum  of  the  along  path 
component  of  the  temperature  gradient  (Mourn  and  Nash,  2009).  We  then  explored  the  use  of  MLE 
spectral  fitting  to  our  data  (Ruddick  et  al.,  2000).  However,  since  our  spectra  do  not  resolve  the  high- 
frequency  part  of  the  spectrum  (up  to  the  Batchelor  cut-off  wavenumber  kB),  the  method  described 
by  Mourn  and  Nash  (2009)  is  more  appropriate  for  our  data.  This  approach  is  based  on  integrating 
the  theoretical  spectrum.  For  consistency,  the  data  are  again  filtered  for  platform  motion,  which  does 
however  not  affect  the  frequency  range  of  the  spectrum  involved  in  the  Batchelor  fit.  Furthermore,  a 
spectral  correction  to  the  signal  measured  by  the  thermistor  is  applied  through  a  single  pole  response 
function 


H2(k)  = - x - , 

[1  +  (Uk/fc)2] 


(5) 


with  fc  =  20  Hz  (Mourn  and  Nash,  2009;  Bogucki  et  al.,  2012).  The  correction  is  done  by  multiplying 
the  spectra  by  the  inverse  of  the  filter  function. 

To  ensure  that  U  varies  as  little  as  possible  over  the  segment  used  to  calculate  temperature  gradient 
spectra,  spectral  estimates  were  done  from  4-s  data  segments  (128  data  points)  and  the  results  for  x 
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Fig.  6.  Flowchart  visualizing  the  processing  steps  used  for  calculation  of  TKED  (top)  and  TD  (bottom). 


were  averaged.  The  length  of  4  s  for  the  data  segments  was  chosen  as  a  trade-off  between  having 
an  adequate  number  of  points  to  calculate  the  spectra,  while  keeping  the  segment  length  as  short 
as  possible  to  limit  variability  of  the  mean  flow.  Our  sampling  frequency  of  32  Hz  is  at  least  three  to 
four  times  lower  than  that  generally  used  to  estimate  /  from  temperature  measurements.  Despite  this 
shortcoming,  which  is  due  in  part  to  the  limitation  of  the  temperature  probe  being  integrated  with  and 
controlled  by  the  ADV,  we  find  the  TD  rates  we  obtain  from  either  method  to  be  robust.  The  calculated 
values  of  TD  are  consistent  between  data  from  different  sensors  at  the  same  measurement  time  and 
depth,  nonetheless,  uncertainties  are  introduced  in  a  number  of  ways  (frame  movement,  choice  of  U, 
mean  temperature  gradient,  among  others),  which  cannot  adequately  be  quantified  at  this  point. 

The  Batchelor  spectral  fitting  technique  can  also  be  applied  to  estimate  e,  through  the  relationship 


which  ties  the  Batchelor  wavenumber  kB  to  the  TKED  rate  e  (v  is  the  molecular  viscosity  of  water),  as 
well  as  the  relationship 


= 


N2X 

2^<s  r 


(7) 


where  ex  signifies  TKED  derived  from  /.  Here,  r  =  0.2,  N  is  the  local  buoyancy  frequency,  and  (^),  is 
the  local  temperature  gradient,  which  can  be  estimated  directly  from  our  data  following  the  procedure 
described  in  Mourn  and  Nash  (2009).  The  mean  temperature  gradient  (^)  is  estimated  from  1-min 
long  data  segments.  Likely  due  to  the  fact  that  our  spectra  are  severely  under  resolved  in  the  range 
preceding  kB,  as  well  as  due  to  uncertainties  associated  with  local  estimates  of  (^)  and  N2,  we  found 
our  estimates  of  ex  from  the  temperature  gradient  spectra  to  be  less  robust  than  those  of  /  itself. 
The  individual  processing  steps  for  TKED  and  TD  are  visualized  in  the  flow  chart  shown  in  Fig.  6. 


2.3.  Optical  measurements  on  the  iMast 

The  velocity  and  temperature  data  on  the  iMast  were  collected  in  the  context  of  optical  target 
clarity  and  spatial  resolution.  A  resolution  chart  was  mounted  to  one  side  of  the  iMast,  a  high- 
definition  camera  to  the  other.  Video  of  the  optical  chart  were  taken  during  the  deployments  and  the 
images  assessed  for  image  degradation.  The  optical  target  used  in  this  analysis  was  the  image  of  an 
optical  chart  on  an  iPad  screen,  hence  an  active  target.  Only  one  pattern  was  used  in  this  study  (of  the 
1951  USAF  resolution  test  chart  type),  but  the  setup  with  the  iPad  as  active  target  allows  the  use  of  a 
variety  of  optical  charts,  including  slant  angle  and  grid  patterns.  The  camera  was  a  Casio  EXILIM  EX-FI, 
which  saves  video  in  MOV  format,  H.264/AVC  standard.  The  videos  were  shot  in  different  resolutions; 
near  the  Florida  coast,  the  video  was  shot  in  High-Speed  (HS)  with  512  by  384  pixels  at  300  frames 
per  second  (fps),  while  in  the  Bahamas,  the  video  setting  was  High-Definition  (HD),  with  60  fps  and 
resolution  1920  by  1080  pixels.  Since  the  data  collection  was  plagued  by  noise  issues  due  to  platform 
motion,  as  described  above,  these  were  the  only  two  videos  where  we  also  had  coinciding  and  usable 
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turbulence  measurements.  For  subsequent  image  analysis,  frames  were  extracted  as  TIFF  images  from 
the  MOV  videos.  To  quantify  the  extent  of  image  degradation  from  optical  turbulence,  a  common 
approach  involving  the  Structural  Similarity  Index  (SSIM)  Method  (Wang  et  al.,  2004)  was  used.  The 
SSIM  method  is  an  image  quality  metric  and  measures  the  similarity  between  two  images,  when  one 
is  considered  the  “perfect  quality”  reference  image. 

The  SSIM  value  is  calculated  as  the  measure  between  two  windows  x  and  y  of  common  size  N  by 
N  as 


SSIM(x ,  y) 


(2flxfly  +  C\)(2axy  +  C2) 

Ox  +  My  +  cOOx2  +  ay  +  C2)  ' 


(8) 


Here,  /xx  is  the  average  of  x,  pty  is  the  average  ofy,  a2  is  the  variance  of  x,  and  a2  the  variance  ofy. 
<7xy  is  the  covariance  of  x  and  y.  C\  =  (k\l)2  and  c2  =  (k2l)2  are  variables,  where  L  is  the  dynamic 
range  of  the  pixel  values  and  default  values  are  for  k\  =  0.01  and  k2  =  0.03.  The  SSIM  index  then 
becomes  a  value  between  —1  and  1,  where  1  is  the  value  achieved  for  identical  sets  of  data.  A  more 
detailed  description  of  the  method  and  comparison  to  other  image  metrics  such  as  mean  square  error 
(MSE)  and  peak  signal-to-noise  ratio  (PSNR)  can  be  found  in  Wang  et  al.  (2004).  For  our  calculation, we 
used  the  Matlab  script  provided  at  “https://ece.uwaterloo.ca/~z70wang/research/ssim/ 

To  account  for  frame  motion  and  vibration  in  the  videos,  successive  images  were  image  registered 
using  “imregister”  in  Matlab’s  Image  Processing  Toolbox,  which  is  an  intensity-based  image  registra¬ 
tion  algorithm  (see  “http://www.mathworks.com/help/images/ref/imregister.html”).  In  addition  to 
the  SSIM,  we  also  calculated  the  sum  of  absolute  differences  (SAD)  for  subsequent  video  frames.  This 
simple  method  eliminates  the  need  for  a  reference  image.  The  value  of  SAD  was  then  scaled  to  fit  into 
the  range  0-1  and,  for  comparison  with  SSIM,  we  calculated  1  —  SAD. 


3.  Results 


3.1.  TKED 

We  compare  data  from  two  stations,  Station  1  near  the  Florida  coast,  on  the  shelf  in  an  area 
influenced  by  the  Gulf  Stream,  and  Station  3d  in  the  Bahamas,  in  a  deep  basin  called  the  Tongue 
of  the  Ocean  (Fig.  1).  Repeat  casts  from  the  VMP  near  Station  1  reveal  an  energetic  environment 
in  the  upper  ocean,  with  TKED  rates  e  reaching  as  high  as  10-5  m2  s-3  in  the  top  20  m  above  the 
thermocline  and  around  10-7  to  10-6  m2  s~3  in  the  thermocline  (Fig.  7).  The  station  was  occupied 
for  two  days,  June  30  and  July  2,  the  VMP  casts  were  taken  on  July  2.  The  differences  in  mixed  layer 
depth  and  stratification  between  the  two  days  can  be  explained  by  the  temporal  and  spatial  variability 
characteristic  of  this  highly  energetic  area  on  the  east  Florida  shelf,  especially  since  the  boat  was 
drifting  during  deployments.  The  variability  shown  during  successive  VMP  casts  is  also  typical  of 
such  a  deployment  strategy,  since  the  boat  will  drift  in  and  out  of  strongly  turbulent  regions,  which 
in  themselves  are  intermittent  and  vary  in  both  space  and  time.  The  TKED  rates  from  the  ADV  on 
the  iMast  fall  well  within  the  range  of  the  VMP  measurements.  The  ADV  data  was  clean  enough  for 
the  estimation  of  c  at  ~20  m  depth  from  all  three  ADVs  on  both  days,  and  the  values  center  around 
c  ~  10-6  m2  s-3.  Excessive  noise  and  phase-wrapping  in  several  probes  at  deeper  depths  led  to 
fragmented  data  that  did  either  not  resolve  the  ISR  or  exhibited  unrealistically  high  values  of  e  due 
to  noisy  velocity  series.  These  data  had  to  be  rejected.  The  data  from  which  e  could  be  calculated 
for  the  deeper  pauses,  at  ^27  m  on  June  30  and  at  ~43  m  on  July  2  show  that  they  fall  within  the 
range  of  values  shown  by  the  VMP.  The  large  margin  of  error  at  around  27  m  can  also  be  attributed  to 
high  levels  of  noise  and  phase-wraps  in  the  velocity  data.  The  data  from  ~10mon  June  30  show  TKED 
values  significantly  lower  that  those  of  the  VMP  on  July  2.  This  is  not  unexpected  given  the  intermittent 
and  spatially  varying  nature  of  turbulence,  and  can  be  explained  by  different  wind  conditions  for  the 
two  days,  affecting  the  mixing  in  the  upper  ocean,  as  well  as  spatial  differences  between  cast  locations 
taken  from  the  drifting  vessel.  Estimates  of  TKED  from  the  microstructure  temperature  data,  also 
mostly  fall  within  the  range  of  e  measured  by  the  VMP  and  the  ADVs  on  the  iMast,  but  as  explained  in 
Section  2.2.2,  they  were  considerably  less  robust  than  the  values  derived  from  the  ADV  measurements. 
Since  the  frame  of  the  iMast  is  tethered  to  the  drifting  vessel,  and  affected  by  this  drift,  as  well  as 
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Fig.  7.  TKED  rates  from  Station  1  on  the  Florida  coast  from  VMP  and  Nortek  ADV  on  the  iMast  (left).  The  station  was  occupied 
on  June  30  and  July  2,  2011.  Stratification  and  c532  from  the  CTD/optics  package  at  this  station  are  shown  on  the  right.  The 
thin  lines  show  different  realizations  of  VMP  drops  for  e  on  July  2.  Estimates  of  TKED  rate  e  from  the  ADV  data  are  shown  as 
stars  (June  30)  and  circles  (July  2)  with  associated  error  bars  (2  standard  errors  estimated  via  MLE  technique).  Individual  points 
signify  data  from  three  different  instruments,  color-coded  red,  purple  and  blue.  Diamonds  (July  2)  and  squares  (June  30)  are 
estimates  of  ex  from  temperature  gradient  spectra  from  the  microstructure  temperature  probe  on  the  iMast. 


cable  motion  due  to  swell  and  vessel  movement,  the  frame  is  expected  to  generate  a  non-negligible 
amount  of  turbulence  while  being  dragged  through  the  water.  Thus,  we  would  generally  expect  the 
values  of  turbulence  parameters  from  the  measurements  on  the  iMast  to  be  higher  than  or  on  the  high 
end  of  those  from  the  VMP,  except  for  the  case  of  very  calm  conditions,  such  as  may  be  encountered 
in  lakes  or  estuaries  (Hou  et  al.,  2012b).  Nevertheless,  it  is  important  to  be  able  to  characterize  the 
turbulence  experienced  by  the  optics  on  the  iMast  itself,  and  not  merely  compare  the  optics  data  to 
the  background  turbulence  data  from  the  VMP.  The  data  we  collected  was  able  to  serve  this  purpose. 

The  environment  in  the  Tongue  of  the  Ocean,  a  stretch  of  deep  water  in  the  central  Bahamas, 
near  the  Berry  Islands,  is  found  to  be  considerably  more  quiescent  than  the  energetic  ambient  on 
the  east  Florida  shelf.  TKED  rates  found  with  the  VMP  in  this  area  on  July  7  center  around  e  ~ 
10-9-10-8  m2  s-3,  with  values  reaching  as  high  as  10-7  m2  s~3  in  the  thermocline  and  10-6  m2  s-3 
closer  to  the  ocean  surface  (Fig.  8).  The  values  from  the  ADVs,  where  all  pauses  were  taken  around 
50  m  depth,  are  higher  than  those  of  the  VMP,  and  center  around  10~6  m2  s-3.  The  estimates  of  ex 
from  the  microstructure  temperature  data  also  fall  within  this  range.  The  data  at  this  station  was  very 
noisy  and  phase-wrapping  was  strong,  to  the  point  where  the  ADV  data  from  the  repeat  occupations 
of  this  station  on  July  6  and  July  8  were  unusable.  Video  frames  from  the  camera  on  the  iMast  taken 
during  deployments  revealed  vibrations  of  the  frame  and/or  camera  housing.  These  vibrations  may 
be  the  cause  of  the  strong  noise  and  excessive  phase-wrapping  found  in  the  data  from  this  station. 
This  may  also  have  affected  the  data  on  July  7  from  which  e  was  calculated,  and  may  have  led  to  a 
high  bias  in  the  results.  An  alternative  explanation  for  the  high  values  from  the  ADV/CTs  at  this  depth 
and  station,  is  that  the  frame  on  which  the  probes  were  mounted  caused  an  increase  in  ambient  tur¬ 
bulence  levels  in  this  quiescent  environment.  The  frame,  as  it  was  pulled  through  the  water  by  the 
drifting  boat,  reached  vertical  velocities  up  to  ~30  cm  s_1  and  variations  in  depth  of  ~30  cm  (Fig.  3). 
The  effect  of  the  frame  on  ambient  turbulence  levels  may  have  been  less  noticeable  in  the  more  en¬ 
ergetic  waters  on  the  Florida  shelf.  Likely  the  strongest  impact  was  due  to  the  ambient  atmospheric 
and  oceanographic  conditions;  weather  conditions  on  the  Florida  shelf  were  calm,  yet  with  strong 
drift  currents,  whereas  conditions  during  the  Tongue  of  the  Ocean  deployments  were  breezy  with 
significant  swell,  which  may  have  contributed  to  cable  and  frame  motion. 


3.2.  TD 

Estimates  of  TD  rates  x  from  the  microstructure  temperature  probe  on  the  iMast  exhibit  a  similar 
pattern.  The  values  largely  fall  within  the  range  predicted  by  the  VMP,  except  at  the  Bahamas  station, 
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Fig.  8.  TKED  rates  from  Station  3  (Bahamas,  Tongue  of  the  Ocean)  from  VMP  and  Nortek  ADV  on  the  iMast  (left).  Stratification 
and  c532  from  the  CTD/optics  package  at  this  station  are  shown  on  the  right.  The  thin  lines  show  different  realizations  of 
VMP  drops  for  e.  Estimates  of  TKED  rate  e  from  the  ADV  data  are  shown  as  circles  with  associated  error  bars  (2  standard  errors 
estimated  via  MLE  technique).  Individual  points  signify  data  from  three  different  instruments,  color-coded  red,  purple  and  blue. 
Diamonds  are  estimates  of  ex  from  temperature  gradient  spectra  from  the  microstructure  temperature  probe  on  the  iMast. 


Fig.  9.  TD  rates  from  Station  1  (FL  coast)  and  Station  3  (Bahamas,  Tongue  of  the  Ocean)  from  VMP  and  temperature 
microstructure  probe  on  the  iMast.  The  thin  lines  show  different  realizations  of  the  VMP  drops,  the  diamonds  are  TD  estimates 
from  temperature  gradient  spectra  from  the  microstructure  temperature  probe.  Squares  are  for  a  repeat  realization  of  Station 
SI  (from  June  30, 2011).  The  station  was  occupied  on  June  30  and  July  2, 201 1.  Individual  points  signify  data  from  two  different 
instruments,  color-coded  red  and  purple. 


and  for  the  pause  above  the  thermocline  at  the  Florida  station  (Fig.  9).  Especially  at  the  Florida  station, 
the  TD  rates  from  the  VMP  can  be  seen  to  show  a  significant  increase  in  the  thermocline,  up  to  three 
orders  of  magnitude  larger  than  in  the  well-mixed  unstratified  near-surface  layer.  The  Tongue  of  the 
Ocean  stations  also  exhibits  an  increase  in  TD  near  the  thermocline,  but  it  is  less  pronounced  than  at 
the  Florida  station,  consistent  with  a  much  weaker  thermocline  near  Station  3d.  The  values  of  /  from 
the  CT  data  near  20  m  depth  in  the  Florida  stations,  which  are  slightly  higher  than  those  of  the  VMP, 
can  be  explained  by  the  vessel  drifting  northwestward  and  onto  the  shelf  into  shallower  water  during 
deployment.  The  iMast  data  at  Station  1  was  taken  after  the  VMP  casts,  and  the  drift  direction  of  the 
vessel  onto  the  narrow  shelf  would  have  led  the  probes  to  sample  through  an  area  right  above  the 
shoaling  thermocline,  showing  increased  temperature  dissipation. 

The  estimates  of  x  in  the  Bahamas  from  the  iMast  data  at  around  50  m  depth  also  show  values 
near  or  above  the  upper  bound  of  those  seen  by  the  VMP.  The  time  series  were  taken  right  above 
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the  thermocline  in  this  area,  and  here,  swell-induced  vertical  motion  of  the  frame  across  layers  of 
increased  TD,  or  even  mixing  induced  by  the  frame  itself,  could  certainly  have  contributed  to  the 
slight  high  bias  in  the  /  estimates.  However,  instrument  vibration  and  noise  contamination  cannot 
be  completely  ruled  out  and  may  also  have  affected  the  estimates  of  x  at  this  station,  especially 
since  high-frequency  noise  was  visible  in  these  temperature  time  series.  Nonetheless,  note  that  the 
X  estimates  from  the  microstructure  temperature  data  in  the  Tongue  of  the  Ocean,  which  were  taken 
just  above  the  thermocline,  center  around  10-7  C2  s-1 ,  about  an  order  of  magnitude  smaller  than  those 
above  or  near  the  thermocline  at  the  Florida  station,  consistent  with  the  difference  in  turbulence  levels 
seen  at  the  two  stations. 

3.3.  Effect  on  optics-Sn 

Since  the  main  focus  of  this  paper  is  to  describe  the  turbulence  parameters  collected  with  the  in¬ 
struments  on  the  iMast,  we  only  briefly  describe  the  optical  measurements,  which  were  taken  to  as¬ 
sess  the  effect  of  turbulence  on  image  degradation.  The  effect  of  turbulence  on  the  optical  transfer 
function  can  be  described  as  a  function  of  TKED  and  TD,  since  the  optical  turbulence  coefficient  can 

be  defined  as  Sn  ~  (Hou,  2009). 

Sn  is  a  function  of  microstructure  spatial  frequency  and  imaging  range  and  enters  the  equation  for 
the  optical  transfer  function  (OTF)  for  the  case  of  underwater  imaging,  which  has  been  described  as 
follows  (Hou,  2009),  accounting  for  path  radiance,  particle,  and  turbulence  scattering 


OTF (ty ,  r) total  —  OTF (*F ,  v)path  •  OTF (*F ,  v)particle  ‘  OTF (*F ,  r) turbulence 


(9) 


where  is  the  angular  spatial  frequency,  r  is  the  imaging  range,  D  is  the  normalized  radiance,  b  is  the 
scattering  coefficient,  00  relates  to  the  mean  scattering  angle,  and  c  is  the  beam  attenuation  coefficient 
for  a  desired  wavelength.  Here,  we  adopt  this  for  the  broad  spectrum,  as  it  only  affects  high-order 
terms  (see  Hou  et  al.,  2007). 

Thus,  Sn  is  directly  related  to  and  has  been  derived  from  the  parameter  commonly  used  in  the 
investigation  of  atmospheric  optical  turbulence,  C2,  see  Hou  (2009).  The  units  for  this  parameter  in 
the  literature  are  based  on  semi-empirical  relations  and  have  no  clear  physical  meaning  (Roggemann 
and  Welsh,  1996).  Sn  is  given  with  the  units  of  rad-  5  m-1,  when  it  appears  in  the  equation  for  OTF, 
however,  when  calculated  as  Sn  ~  xe-3 ,  the  units  are  neglected  (Hou  et  al.,  2012b).  More  details  on 
Sn  and  its  relationship  to  the  optical  transfer  function  can  be  found  in  Hou  (2009);  Hou  et  al.  (2012b). 

We  estimate  values  of  Sn  for  Stations  1  and  3d,  in  Florida  and  the  Bahamas,  respectively,  for  the 
pauses  near  or  just  above  the  thermocline,  where  we  can  expect  the  strongest  optical  turbulence. 

At  Station  1  near  the  Florida  coast,  at  ~20  m  depth,  with  a  beam  attenuation  at  532  nm  of 
c532  «  0.15  m-1,  using  532  nm  as  a  proxy  (Hou  et  al.,  2007),  the  value  becomes 

Sn  -  10-7  C2  s-1(10-6m2  s-3)"5  ^  10-5  with  the  VMP  data,  or 

Sn  -  10-6  C2  s-1(10-6  m2  s-3)"5  «  10-4  from  ADV/CT  data  on  the  iMast  (Table  1). 

In  the  Bahamas,  near  50  m  depth  and  with  a  c532  ~  0.1  m-1, 

Sn  ~  10-9  C2  s-1(10-8  m2  s-3)"5  «  5  X  10-7  from  the  VMP  and 

Sn  ~  10-7  C2  s-1(10-6  m2  s-3)“3  «  10-5  from  the  iMast. 

The  data  from  the  iMast  show  higher  values  of  x  and  e  for  both  stations,  which  is  not  unexpected, 
given  the  complexities  and  challenges  involved  in  interpreting  the  ADV  data  collected  on  a  platform 
subject  to  wave  motion  and  cable  vibration  (Table  1).  Nonetheless,  the  data  are  able  to  capture  the 
trend. 

Not  surprisingly,  the  impact  of  optical  turbulence  is  found  to  be  more  pronounced  in  the  more 
energetic  environment  near  the  Florida  coast  and  less  apparent  in  the  deeper,  quiescent  waters  in 
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Table  1 

Optical  turbulence  parameter  S„. 


Location 

Instrument 

Depth  (m)  c532(m 

X  (C2  s 

')  €(m2s-3) 

Sn 

FL 

VMP 

20 

0.15 

10~7 

10“6 

10“5 

FL 

iMast 

20 

0.15 

10“6 

10“6 

10“4 

Bah 

VMP 

50 

0.1 

10-9 

10“8 

5 ■ 10-7 

Bah 

iMast 

50 

0.1 

10-7 

10“6 

10“5 

Table  2 

Optical  transfer  function  OTF. 

Location 

Instrument 

c532  (m-1) 

Sn 

OTF  particle 

OTF  turbulence 

FL 

VMP 

0.15 

10-5 

0.86 

0.37 

FL 

iMast 

0.15 

10~4 

0.86 

0.0001 

Bah 

VMP 

0.10 

10-6 

0.91 

0.90 

Bah 

iMast 

0.10 

10~5 

0.91 

0.37 

the  Tongue  of  the  Ocean.  Frames  taken  from  the  videos  on  the  iMast  confirm  this  finding  (Figs.  10- 
13).  Frames  are  shown  at  10  s  intervals,  and  illustrate  the  variability  over  the  course  of  about  1  min, 
especially  near  the  Florida  coast.  At  the  same  level  of  magnification,  the  image  from  the  Bahamas 
stations  appears  clearer,  especially  at  high  frequencies,  than  the  image  from  the  Florida  station.  Near 
the  Florida  Coast,  the  variability  is  particularly  pronounced  at  higher  spatial  frequencies,  indicative  of 
the  effects  of  turbulent  density  variations  on  the  images.  Note  that  beam  c  enters  the  equation  for  the 
optical  transfer  function  (Eq.  (9))  in  the  exponent,  thus,  while  the  beam  c  values  measured  at  the  two 
stations,  indicative  of  levels  of  particle  concentrations,  differ  slightly,  the  0.05  m-1  difference  in  beam 
c  cannot  account  for  all  of  this  image  degradation.  This  can  be  seen  from  a  simple  calculation  of  OTFs 
using  Eq.  (9).  We  can  choose  00  =  0.03,  r  =  1  and  V  ~  1000,  as  well  as  b  =  0.12  m-1  for  Florida 
and  b  =  0.08  m-1  for  the  Bahamas.  This  will  result  in  OTFparticie  =  0.86  for  a  c  =  0.15  m-1  as  seen 
at  the  Florida  station,  and  OTFparticie  =  0.91  for  a  c  =  0.10  m-1  as  seen  at  the  Bahamas  station.  The 
contribution  of  the  turbulence,  when  calculated  using  the  values  for  /  and  e  from  the  VMP,  comes 
to  OTF turbulence  =  0.37  for  S„MP  =  10-5,  such  as  at  the  Florida  station,  and  OTF turbulence  =  0.90  for 
svmp  _  jo-6,  for  the  Bahamas  Station.  Table  2  summarizes  these  values.  As  noted  previously,  there 
are  some  uncertainties  associated  with  the  numbers  calculated  from  the  ADV/CT  data  on  the  frame, 
but  the  data  are  able  to  capture  the  general  trend.  It  is  interesting  to  note  that,  when  optical  turbulence 
is  extremely  high,  Sn  ~  10-4,  then  OTFturb  tends  to  zero  for  high  spatial  frequencies,  which  is  in  line 
with  the  Simple  Underwater  Imaging  Model  (SUIM)  described  in  FIou  (2009). 

We  acknowledge  that  the  particle  fields  at  the  two  locations  are  also  likely  to  differ  and  thus  affect 
the  images,  however,  the  detailed  investigation  of  this  influence  is  beyond  the  scope  of  this  paper. 

To  quantify  the  extent  of  image  degradation,  including  imaged  degradation  due  to  the  effect  of 
turbulence  microstructure,  we  apply  the  Structural  Similarity  Index  Method  (SSIM)  (Wang  et  al.,  2004) 
to  a  subset  of  the  images  containing  the  higher  spatial  frequencies,  since  the  effect  of  turbulence  on 
image  degradation  is  especially  pronounced  at  higher  spatial  frequencies. 

Here,  we  use  images  taken  near  the  water  surface,  shortly  after  the  frame  enters  the  water  and  is 
held  steady  near  the  surface,  as  the  respective  reference  images  for  the  two  stations.  This  also  serves 
to  effectively  normalize  the  particle  fields  and  ambient  light  conditions  for  the  respective  stations. 

The  SSIM  values  are  indicated  in  Figs.  1 1  and  13  and  shown  as  a  time  series  in  Fig.  14.  They  confirm 
that  the  image  degradation  as  well  as  the  low-frequency  temporal  variability,  which  is  indicative  of 
turbulence  effects,  are  greater  near  the  Florida  coast.  The  high-frequency  temporal  variability  and 
larger  spread  seen  in  the  images  from  the  Bahamas  is  related  to  vibrations  of  the  frame,  likely  due  to 
cable  motion  enhanced  by  the  sea-state.  These  vibrations  can  be  seen  in  the  video  images  and  may 
not  be  entirely  removed  by  the  shift  invariance  algorithm.  An  alternative  method  to  the  SSIM,  is  to 
calculate  the  SAD  for  subsequent  image  frames.  The  resulting  time  series  is  also  shown  in  Fig.  14. 
The  SAD  values  (normalized  to  the  range  of  0-1  and  shown  as  1  —  SAD)  overall  follow  the  SSIM, 
however,  there  appears  to  be  increased  high-frequency  variability  in  the  Bahamas  station,  especially 
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Fig.  10.  Frames  from  the  HS  video  taken  on  the  iMast  near  the  Florida  coast  (Station  SI  at  20  m  depth;  resolution  512  by 
384  pixels  at  300  fps)  are  shown  at  10  s  intervals  (time  progression  is  from  left  to  right  and  top  to  bottom,  starting  with  top 
right  image).  Reference  image  is  shown  on  the  top  left.  The  dashed  rectangle  indicates  the  part  of  the  image  that  is  shown  in 
more  detail  in  Fig.  11. 


towards  the  beginning  of  the  time  series.  This  can  be  attributed  to  the  aforementioned  high-frequency 
vibrations  seen  in  the  images  from  this  station.  Any  lateral  shift  in  the  image,  such  as  caused  by 
vibrations  will  lead  to  a  high  value  of  SAD,  whereas  turbulence,  which  leads  to  the  deformation  of 
adjacent  pixels,  will  be  better  detected  by  the  SSIM  algorithm. 
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Fig.  11.  Close-up  of  the  area  of  the  dashed  rectangle  from  Fig.  1 0  highlights  the  variability  at  higher  spatial  frequencies,  at  10  s 
intervals  (time  progression  is  from  left  to  right  and  top  to  bottom,  starting  with  image  on  top  right).  The  black  numbers  are 
the  values  from  applying  the  SSIM  algorithm  to  these  images  using  a  reference  image  (top  left)  from  near  the  surface  at  this 
location. 


4.  Summary  and  discussion 

In  order  to  test  a  new  platform  to  study  the  impact  of  oceanic  optical  turbulence  on  EO  signal 
transmission,  we  collected  high-resolution  velocity  and  temperature  data  along  with  high-resolution 
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Fig.  12.  Frames  from  the  HD  video  taken  on  the  iMast  in  the  Bahamas  (Station  S3  at  50  m  depth;  resolution  1920  by  1080  pixels 
at  60  fps)  are  shown  at  10  s  intervals  (time  progression  is  from  left  to  right  and  top  to  bottom,  starting  with  top  right  image). 
Reference  image  is  shown  on  the  top  left.  The  dashed  rectangle  indicates  the  part  of  the  image  that  is  shown  in  more  detail  in 

Fig.  13. 


video  frames  in  the  waters  off  the  Florida  coast  and  the  Bahamas.  A  VMP  was  deployed  to  collect 
profiles  of  velocity  shear  and  temperature  data.  TKED  and  TD  rates  were  inferred  from  the  ADV  and 
temperature  microstructure  data,  compared  to  the  dissipation  rates  from  the  VMP  and  put  into  the 
context  of  the  optical  measurements. 

There  are  considerable  inherent  challenges  to  the  study  of  optical  turbulence  in  the  open 
ocean.  The  temporal  and  spatial  scales  are  vastly  different  between  the  quasi-instantaneous  optical 
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Fig.  13.  Close-up  of  the  area  of  the  dashed  rectangle  from  Fig.  12,  at  10  s  intervals  (time  progression  is  from  left  to  right  and 
top  to  bottom,  starting  with  image  on  top  right).  The  black  numbers  are  the  values  from  applying  the  SSIM  algorithm  to  these 
images  (using  a  reference  image  from  near  the  surface  at  this  location,  shown  on  top  left)  and  show  less  variability  than  in 

Fig.  11. 

measurements  (on  the  order  of  ms)  and  the  statistical  averages  needed  for  velocity  and  temperature 
data  to  estimate  turbulence  properties  (on  the  order  of  tens  of  minutes).  This  scale  discrepancy  is 
compounded  by  the  temporal  and  spatial  variability  of  oceanic  turbulence,  which  is  well-known  to  be 
intermittent  and  “patchy”. 

Measuring  turbulence  in  the  ocean  is  a  challenging  task,  complicated  by  deployment  strategy,  ap¬ 
parent  in  our  data  as  the  drift  across  a  shoaling  thermocline  on  the  Florida  shelf,  and  likely  frame- 


56 


S.  Matt  et  al.  /  Methods  in  Oceanography  1  1  (2014)  39-58 
Florida  Station,  SI 


0.2  - 
0.1 

O' - 1 - *- 

0  0.5  1 


1,5  2  25 

time  (min) 


Bahamas  Station,  S3 


SSIM  w/  ref  img 

- SSIM  running  mean 

- SAD  running  mean,  no  ref  img 


0  0.5  1  1.5  2  2.5 

time  (min) 

Fig.  14.  Time  series  of  the  SSIM  algorithm  applied  to  video  stills  from  near  the  Florida  Coast  (top)  and  the  Bahamas  (bottom) 
over  several  minutes  illustrate  the  variability  in  image  degradation  at  both  locations.  The  gray  line  shows  the  running  mean  of 
the  two  time  series.  The  low-frequency  temporal  variability  and  the  image  degradation  are  more  pronounced  near  the  Florida 
Coast.  The  dashed  red  line  shows  the  running  mean  of  the  sum  of  absolute  differences  (SAD)  for  subsequent  video  frames  in 
the  time  series.  The  value  of  SAD  is  scaled  to  fit  into  the  range  0-1  and  shown  as  1  —  SAD.  (For  interpretation  of  the  references 
to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


induced  turbulence  in  the  Bahamas,  as  well  as  the  ambient  environmental  conditions  of  ocean  and 
atmosphere.  Furthermore,  the  use  of  ADVs  to  infer  TKED  is  inherently  more  successful  in  more  ener¬ 
getic  environments,  such  as  that  encountered  on  the  Florida  shelf.  In  the  calmer  waters  of  the  Tongue 
of  the  Ocean,  Bahamas,  the  signal  to  noise  ratio  of  the  sensor  may  get  close  to  its  lower  limit.  The  noise 
from  platform  motion  and  vibrations,  apparent  mainly  in  our  velocity  data,  can  further  complicate  the 
estimation  of  turbulence  properties  from  the  ADV/CT  package.  Thus,  complementary  measurements 
of  ambient  turbulence  levels  with  a  dedicated  microstructure  profiler,  such  as  the  VMP,  are  an  invalu¬ 
able  tool  to  aid  in  data  analysis  and  guide  error  estimates. 

Despite  these  difficulties,  the  data  from  ADV  and  temperature  microstructure  probe  were  able  to 
provide  information  on  the  turbulence  experienced  by  the  iMast,  while  simultaneously  collecting  in¬ 
formation  on  the  impact  of  this  turbulence  on  image  degradation.  Increased  levels  of  optical  turbu- 
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lence  were  expected  to  be  found  above  or  near  the  thermocline  and  these  optical  turbulence  levels 
predicted  by  the  VMP  and  ADV/CT  data  were  confirmed  by  the  images  taken  during  the  iMast  deploy¬ 
ments.  This  study  was  the  first  comprehensive  effort  collecting  data  from  both  a  traditional  profiling 
approach  for  measuring  turbulence,  as  well  from  an  ADV/CT  package  on  a  moving  platform,  for  the 
investigation  of  optical  turbulence  in  the  energetic  open  ocean  environment. 

While  there  are  many  challenges,  the  effort  of  studying  optical  turbulence  in  the  natural  environ¬ 
ment  is  necessary  to  quantify  the  impact  of  oceanic  microstructure  on  underwater  EO  signal  trans¬ 
mission  and  to  advance  our  knowledge  on  the  use  of  optical  techniques  to  measure  turbulence.  Future 
field  studies  of  optical  turbulence  may  benefit  from  a  platform  design  that  allows  the  frame  to  freely 
drift  or  fall  through  the  water  column,  while  still  maintaining  an  adequate  optical  path  length. 

The  results  of  our  study  will  help  address  the  sampling  challenges  for  future  field  efforts  aimed  at 
quantifying  optical  turbulence  and  may  also  improve  our  understanding  of  the  use  of  ADVs  to  probe 
oceanic  turbulence.  In  addition  to  studies  of  naturally-occurring  optical  turbulence,  laboratory  exper¬ 
iments,  conducted  in  a  controlled  laboratory  environment  and  complemented  by  numerical  simula¬ 
tions,  should  prove  very  useful  to  the  detailed  and  quantitative  investigation  of  the  processes  involved 
in  optical  turbulence  in  the  ocean. 
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